DESY 09-041 
LPTENS-09/09 
AEI-2009-035 



Exact Spectrum of Anomalous Dimensions of Planar = 4 
Supersymmetric Yang-Mills Theory: TBA and excited states 

Nikolay Gromov'^ Vladimir Kazakov* Andrii Kozak° Pedro Vieira® 

'^DESY Theory, Hamburg, Germany & 
II. Institut fiir Theoretische Physik Universitat, Hamburg, Germany & 
St. Petersburg INF, St. Petersburg, Russia 

•'°Ecole Normale Superieure, LPT, 75231 Paris CEDEX-5, France & 
I'Universite Paris- VI, Paris, France; 

®Max-Planck-Institut fiir Gravitationphysik 
Albert-Einstein-Institut, 14476 Potsdam, Germany 



Abstract 

Using the thermodynamic Bethe ansatz method we derive an infinite set of integral non-linear 
equations for the spectrum of states/operators in AdS/CFT. The Y-system conjectured in [1] 
for the spectrum of all operators in planar iV = 4 SYM theory follows from these equations. In 
particular, we present the integral TBA type equations for the spectrum of all operators within 
the sl{2) sector. We prove that all the kernels and free terms entering these TBA equations are 
real and have nice fusion properties in the relevant mirror kinematics. We find the analogue of 
DHM formula for the dressing kernel in the mirror kinematics. 
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Figure 1: T-shaped "fat hook" (T-hook) uniting two SU{2\2) fat hooks, see [2] for details on fat 
hooks and super algebras. 

1 Introduction 



Recently, a set of functional equations, the so called y-system, defining the spectrum of all local 
operators in planar AdS/CFT correspondence, was proposed by three of the current authors [i]. 
The Y-system has the form of functional equations 



^ a,s ^ a,s 



[l + Ya,s+i){l + Ya,s-i) 



a+l,s a—l,s 



(1) 



{l + Ya+l,s)il + Ya-l,s) 

where = f{u±i/2) are simple shifts in the imaginary direction. The functions are 
defined only on the nodes marked by o, e, ®, £\, • on Fig.l. Its solutions with appropriate analytical 
properties define the energy of a state (anomalous dimension of an operator in N=4 SYM) through 
the formula^ 

E^iK.)+f; r ^ w)- (2) 



E 



where e* is the mirror "momentum" defined in the text below and the rapidities u^j are fixed by 
the exact Bethe ansatz equations 

Yi,o{u4,j) = -I ■ (3) 



The Y-system is equivalent to the Hirota bilinear equation 

7"* -1- 



(4) 



where the functions Ta^s{u) are non-zero only on the visible part of the 2D lattice drawn on Fig.l 
and 

^a,s+l^a,s— 1 /r\ 

^ 0+1,5-^0—1,3 

It was shown that the Y-system passes a few non-trivial tests, and in particular it is completely 
consistent with the asymptotic Bethe ansatz (ABA) [I, 5, 6, 7], is compatible with the crossing 



^In some cases the integration contour could encircle singularities of the integrand situated away from the real 
axe. In the large L asymptotics these singularities can be responsible for the Liischer ^-terms. See also discussion in 
section 7. 
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relation [3] and reproduces the first wrapping corrections at weak coupling for Konishi and other 
twist two operators [S, 9, 10]. 

In this paper, we will provide a derivation of the y-system similar in spirit to that employed 
in the derivation of the TBA-type non-linear integral equations for the finite volume spectra of 
relativistic 2-dimensional models. It is based on the Matsubara trick relating the ground state of 
a euclidean QFT on a cylinder to the free energy of the same theory in finite temperature. If we 
take instead of the cylinder a torus with a small circumference L and a large circumference R we 
can represent the partition function in two different channels as a sum over energiy levels. In the 
large R limit, we can identify the free energy J-{L) per unit length of a "mirror" QFT living in the 
space section along the infinite direction of the torus and having a temperature T = 1/L, with the 
ground state energy Eq{L) of the original QFT living on a space circle of the radius L 

Z{L,R) = ^e-^^^(«) =Y^^-RE,iL) ^-RTiL) = ^-RE,(L) 

k j 



In the relativistic QFT's the original theory and the mirror theory are essentially equivalent and 
differ only in the boundary conditions [ ]. An example of such a TBA calculation, useful for our 
further purposes, for the SU{2) principle chiral field (PCF), can be seen in the Appendix A of [ -]. 
In the supersting sigma model on AdS^ x background in the light cone gauge relevant to our 
problem, we have to deal with the non-relativistic original and mirror sigma models (see [1-3, 14]). 

Particularly important for our discussion is the form of the energy and momentum of the 
elementary excitations for both the physical and mirror theories in infinite volume. They are 
conveniently parameterized in terms of the Zhukowsky variables, 

1 u 

x{u) g 

which admits two solutions, one of them outside the unit circle |x(n)| > 1 and another inside the 
unit circle, |x(ii)| < 1. The energy Ca ("u) and momentum ("u) of the physical bound states are 
then given by [ ] 

2ig , X 1^ x^^"^ 
^«(^) = « + ^ - ^ ' P^(^) = z ^ 

where x'^"! = x{u ± ia/2) are evaluated in the physical kinematics where > 1. 

The mirror energy and momentum are obtained by the usual Wick rotation {E,p) {ip,iE). 
To stress this we denote the mirror energy by ip* and the mirror momentum by ie*. The quantities 
e* and are defined precisely as in (7) where x^^^ are now evaluated in the mirror kinematics where 
> 1 but < 1, for a > 0. 

Let us now return to our general review of the TBA method. This method is based on the 
so called string hypothesis: all the eigenstates of an integrable model in the infinite volume are 
represented by bound states (the simplest ones are called "strings") described by some density pA- 
In terms of these densities the asymptotic Bethe equations simply read 

Pa{u) + pa{u) = J_^!aM _ Kba{v,u) * pb{v) . (8) 
zvr au 

Here Kba{viU) = -^-^log Sab{u^v) is the kernel describing the interaction between the bound 
states A and B which scatter via an S-matrix Sab- i(*A is the momentum of a magnon labeled by 
A. For the same reasons as mentioned above in the discussion of the AdS/CFT dispersion relations 
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we use this notation to emphasize that the momenta of these mirror particles are obtained from 
the energy of the physical particles CAiu) by the Wick rotation. Finally is the density of holes 
associated with the bound state A. 



To compute the free energy we must minimize the functional 

-F= Y.j°^ du(^{Lip\ + hA)pA- /?Alog(l + ^) +PAlog(l + ^) ) 



(9) 

with respect to pa{u)-, Pa{u) and exclude 5pA by the use of the constraint imposed by the BAE's (8). 
The physical origin of each term in the expression for the free energy is as follows: The first term 
accounts for the energy (times inverse "temperature" L); the term in the square brackets represent 
the entropy contribution; we added a generic chemical potential Ka for each kind of bound states. 
This chemical potential is needed if the theory contains fermionic excitations, as is the case for the 
AdS/CFT system, since we want to compute the Witten index rather than the thermal partition 
function where the physical fermions are periodic. This amounts to choosing /i^i = ivr = log(— 1) 
for the fermionic states and /i^ = for the bosonic states. 

The minimization of the free energy yields the TBA equations 

\ogyA{u) = Kab{u,v) * log[l + l/yB{v)] + iLp\{u) + hA (10) 

for the quantities yA = Finally, at this saddle point, the free energy can be simply written as 

In this way one obtains the finite volume ground state energy for a generic integrable field theory. 
The excited physical states are recovered by the usual procedure of analytic continuation [16, 17, 
18, 10] and will be also discussed in this paper. 

In what follows, we will apply the TBA method to the "mirror" superstring sigma model and 
derive this AdS/CFT Y-system conjectured in [ ]. The actual TBA equations arising as an in- 
termediate step towards the Y-system, may be very useful for the numerical calculations of the 
energies of low-lying states. 



2 The starting point: Beisert-Staudacher equations 



The basis of our derivation of TBA for AdS/CFT are the Beisert-Staudacher (BS) ABA equations 
of [4, o, ~] in their mirror form [20, 14]. We write them in our compact notations, introducing three 
types of Baxter functions 



b[^\u) 



i=i 



Ki _±_ 
-j-j- x{u) 



X 



Qi{u) 



Ki 

W{u-uij) 



[-g)'^^Ri{u)Bi{u). 

(12) 

The index / takes the values I = IL, 2L, 3L or / = IR, 2R, 3R parametrizing the rapidities of the 
left and right SU{2\2) wings of the model, correspondingly. R^^^ and B^^^ with no subscript 
/ correspond to the roots x^j of the middle node and Ri,Bi without supercript (+) or (— ) are 
defined as in (12) with x^ replaced by xj. In these notations the left wings ABA's read: 



Q-rBW 



Q2L QilQsl 



Q2L QilQsl 



U2L,k 



Q-rRM 



(13) 
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with a similar set of equations for the right wing replacing L ^ R. The Bethe equation for the 
middle node for the full AdS/CFT ABA of [~] fix the positions of the n4j- roots from^ 



Q4 



B^lRsl B^rRsrJ 



n 



X 



2 



4,i 



(14) 



M=«4.fe 



for the sl{2) favored grading. The dressing factor is S{u) = Ylj cr{x{u), x^j) where a is the BES 
dressing kernel [i] (see [21] for a nice integral representation of the dressing kernel). 



3 Bound states and TBA equations for the mirror "free energy" 

To write the TBA for the full AdS/CFT, we have to find the BAE's for the densities of all complexes 
of Bethe roots in the infinite volume R = 00. The string hypothesis implies the full description 
of the infinite volume solutions. They are easy to classify: there is only one type of momentum 
carrying complexes, strings in the middle nodes, similar to standard SU{2) strings [15]; the rest 
are the same complexes as found by Takahashi in the Hubbard model [22, 23] (see also [24]). 

As the result, we find that in the large R limit of BAE's the roots regroup into the following 
bound states: 
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? • • 
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L,R 
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= u 












L, R single fermion 


: e± 


L,R 
^3 


= u 












L, R single fermion 


: ®-|- 



where by u we denote the real center of a complex. Thus the index A in formulae (8-11) takes the 
values 

^ = {o±„,®±,®±,fl±„,«„} (15) 
or, in the notation used in [ ] for the points on the T-hook 

A = {(1, ±n), (2, ±2), (1, ±1), (n, ±1), (n, 0)} . (16) 

^This equation is identical to the eq.(6.6) from [14]. The factors of jx" outside of the square brackets can 
be easily reconstructed from the unimodularity of the r.h.s. of (14). We thank the referee for pointing us out this 
misprint which fortunately does not affect any of our results in the previous version of the preprint. 
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Multiplying the Bethe equations along each complex we obtain the fused equations (8) for the 
densities (of particles and holes, pAiu) and pa{u)) of the centers of complexes (10). It is useful to 
introduce the following notation for 



1 1 



(17) 



In particular notice that the Y functions y^.s arrange nicely into a T-shaped form as depicted in 
Fig.l. As shown below, these functions are precisely those appearing in the y-system (1). 

The only complexes which carry energy and momentum are those made out of middle node 
roots 

4 = , Pa = ^A,»„Pn (18) 

where e* and p* are explained after (7). The fused kernels Kab are given by 



Kab 
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Ofn 
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®_|_ 
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+Kn-1 
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+Km-1 


.^(01) 
"im 


®_|_ 


Km—1 








+Km-1 


7? (01) 







-Kn-1 


+Kn-1 


~^Kn—l^m,—l 


7? (01) fo{01) 
l^nm '~^n-2,m 


•n 





"nl 


-r,(10) 


^(10) „(10) 


no ,-0(11) 



where the block entrees of this infinite matrix are defined as 



(19) 



Id, u — V + m 2 

Kn = log T 

2-Ki dv u — V — in 2 



m — 1 n — 1 

2 2 



Knm = ^ ^ K. 



2j+2fc+2 



T— m-1 J n-1 



Snm{n,v) ^ -LAiog^(x±"(n),x±"(^;)) 
Ztti av 



n— 1 m— 1 

2 2 

E E 

_ 1—1 ; — ™.- 



Id, 6(u + m/2 + ij, v - ib/2 + ik) 

log 

^ 21X1 dv b{u — ia/2 + ij, v + ib/2 + ik) 



(20) 
(21) 
(22) 



71—1 m— 1 

2 2 



n-l I m- 



1 d r(u + ia/2 + ij,v — ib/2 + ik) 

log 

2'iri dv r{u — ia/2 + ij, v + ib/2 + i/c) 



(23) 



where 



x(m)-x(t;) l/x(u)-x(f) 
r(n, u) = -= — , b(u,v) = 



(24) 



\/x{v) ' ' ' ' \/x{v) 

In the table above we only wrote the interaction between the complexes of the left SU{2\2) wing, 
between those complexes and the middle node bound states, as well as between the middle node 
bound states themselves. The right wing interaction is of course absolutely identical and the 
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complexes of different wings do not interact. Equations (10) in the notation of (17) then read 
logl®± = +Km^i * log + nf^^ * log(l + Y. J + log(-l) (25) 



zbm 



logye± = -Km^i * log - * log(l + Y.J - log(-l) (26) 

log Ya^^ = * log(l + yQ^,J - * log ^J" (27) 

+ fe+sH„Uiog(i + y, 



iogyo±„ = iv:„-i,™-i*iog(i + i/yo±„) + K„_i*iog— — (28) 

i + i/ r®^ 

iogy.„ = Liog^+ 25„™-7^(l^) + s(li) *iog(i + y.„) (29) 



- 4f * iog(i + i/y® J + * iog(i + Y^j + {n(Z^ + * iog(i + ya^; 

- sif * iog(i + i/ye_ ) + 7^i\°) * iog(i + y« J + (Te^^Q) + ^^l^) * iog(i + y^_^ 



All convolutions are to be understood in the usual sense with the second variable being integrated 
over so that K * f = J dvK(u,v)f{v). Summation over the repeated index m is assumed (m = 
2, . . . , 00 for the convolutions involving pyramids fl±m and strings o±m and m = 1, . . . , oo for the 
convolutions with the middle node bound states •m)- There are still some ambiguities involved in 
these integral equations concerning the choice of the integration contours. We will discuss this, still 
not completely elucidated, point when we will consider equations for the excited states where some 
of the ambiguities will be lifted. 

4 Derivation of the AdS/CFT Y-system 

We will now derive, from the TBA equations, the Y-system (1) and (4) for the AdS/CFT spectrum 
conjectured in [ ]. We shall do it separately for each type of excitations. 

The key idea in the derivation is to use the discrete Laplace operator acting on the free variable 
u and free index n in the TBA equations. We notice that 

AKn{u) = Kn{u + i/2 - iO) + Kn{u - i/2 + «0) - i^„+i(n) - Kn-i{u) = 5n,i5{u) 

As a simple consequence of this identity we find 

AKnm{v -U) = An'^^{v,u) = 6n,m+lS{v - u) + 6n,m-lS{v - u) 

A7^io^)(^;,7/) = Aiz'-^iv^u) = 5n,m5{v - u) (30) 

whereas the Laplacian kills all other kernels, ASnm = 0, etc. For example, the fact that the dressing 
factor is killed by the Laplacian follows from its harmonic form 

O'nmiu v) = QX{-^+i"-/2,v+in/2)+x{u-in/2,v-in/2)-x{u-in/2,v+in/2)-x{u+in/2,v-in/2) ^g-j^-j 
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without any singularities in the physical kinematics (this fact was aheady used in ['] when con- 
structing the large L solutions of the y-system) . By virtue of these identities we can easily compute 

the combinations log y °" °" — , log y J^" — and log y *" *" — , where = f{u ± i/2 ^ iO), 

using respectively (28), (27) and (29). We find 

log = log(l + l/yo„+J(l + l/>o_J ,n>2 (32) 

and 

^°^^ = ^°^ — tttav — 

for the string bound states. The equations for at n < —2, as well as their derivation, are 
similar. For the pyramid complexes we obtain 

Yn Yn 1 + y. 

log y V =l°g n , y vily ) ' (34) 

and 

log ^ = log n VI ^ " ^®+^®+ + Z^(^ni -Ki )* log(l + ^.J • 

The first term in the r.h.s. of this equation reproduces again the correct structure of the Y-system 
(1). In fact, we will see below that the last two terms cancel each other and hence this equation 
perfectly fits the Y-system (1). Finally, for the middle node bound states, we kill again the kernels 
when applying the discrete Laplace operator and obtain 

log *" *" = log , n > 1 35 

^y.„+iy.„_i ^(i + y.„+j(i + y.„_,) ' ^ ^ 

and 

log V = log -j— T7^ • (36) 

We are left with the equations for the two fermionic nodes Yi^i = Y"®^ and 12,2 = (for 11,-1 
and Y2-2 it will be similar). We consider first the node ^1,1. Combining equations (25) for 
u — )■ n ± i/2 =p zO with equations (27) and (28) for real u and n = 2 we obtain (again using the 
fusion properties of several kernels), 

perfectly reproducing the the equation for Yi^i from the y-system (1). Finally, to find the equation 
for the last fermion node ¥2^2 we simply add up equations (26) and (25) to get 

log y^^ye^ = - Sf™ ) * log(l + Y.J (38) 

m 

This shows indeed that the two last terms in (35) cancel. The equation for Y22 = Y^^ is not a 
part of y-system (1) since in the standard form it would contain the ratio j^^'^Ji?^ = It is thus 
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natural that one can not render this equation local if we only use the finite Y functions, see also 
[25]. However, in terms of the T- functions appearing in 5 we believe, and partially checked, that 
Hirota equation 4 is well defined on the full T-shaped fat-hook of figure 1. 

All these equations precisely reproduce the y-system (1) under the identification 

{yo±„,l^e±,^®±,^a±„,^.±„} = {Yi,±n,Y2,±2,Yi^±uYn,±i,Yn,o} (39) 
mentioned in the previous section! 



5 Integral equations for excited states 

In this section we will consider the non- linear integral TBA-type equations for excited states. For 
simplicity we shall consider only the states in the SL{2) sector, corresponding to operators of the 
form tr [D^ Z'^) -\- permutations. Notice that since none of the wings are excited the K-functions 
will have the symmetry Ya^s = Ya-s which also means that Y®^ = Yi^_ = Y^, .... To consider 
such excited states we employ the standard analytic continuation trick [IG, 17, 19] where we pick 
extra singularities in the convolutions with l*^ at the points where l#i(w4,j) = —1. This procedure 
contains some ambiguities and the result should be considered as a conjecture. In this way, the free 
energy (11) becomes (2) while the non-linear integral equations of section (3) are modified by the 
terms in the square brackets 



logy® 

log^e 
^ogYn 



* log + 7^f^) * log(l + Y.J + 



-Km-l * log 



1 + 
l + Y, 



sS°^)*iog(i + y. 



log 
log 



Ri-) 



+ log(-l) 
-log(-l) 



(40) 
(41) 



* log(l + YaJ- Kn-l * log Y^^^ + {^'nrl} + ^i- U) * log(l + Y. 



+ 



2 



n — 3 
2 



^ ^ R(+\u + ik) ^ ^ 



+ 1A 



n-l 
2 



iogyo„ 



logY.^ 



R(~){u + ik) 

— l,m— 1 

*iog(i + i/yo^) + -fc„-i*iog 



B(-){u + ik) 

1 + Y^ 
1 + 1/Y^ 



(42) 
(43) 



Llog + (2Snm - n^i^ + aii^)) * log(l + Y.J + 



n-l 
2 



i^{u + ik) 



(44) 



+ 2 (7^i\°) * iog(i + y®) - * iog(i + 1/ye) + {n^^^ + ^^2) * iog(i + y^j) 



where 



$(u) = - log 52 



BM+R(-y 
B(-)-R(+y 



(45) 



and B and R and S containing the positions of rapidities of the excited states are defined in section 
2. These rapidities are constrained by the exact Bethe equations 



y«i {uaj 



-1 , i = 1,...,M. 



(46) 
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In the convolutions involving the fermionic y- functions 1© and we integrate over v G [—2(7, 2g] 
We found that prescription to be consistent with the asymptotical large L solution of the Y-system 
derived in [J]. In fact as one can see from these integral equations we can think of the two functions 
Yq and 1/1® as two branches of the same function. In this language the convolutions can be 
recasted into some nice S-cycle contour integrals in the x{u) Riemann sheet. This is reminiscent of 
the inversion symmetry in the BS equations which allows one to reduce the seven Bethe equations 
to a smaller set of five equations [ ]. 

An important check of these equations is the limit where L — )• oo. The solution of the Y-system 
in this limit was constructed explicitly in [J]. We checked numerically that for large L our integral 
equations are consistent with the large volume solution. 

6 Physical and mirror choices of branches 

The above system of TBA equations should be valid for any value of the spectral parameter u 
and it should be possible to analytically continue it to any point of the Riemann surface of the 
multi-valued y-functions. But the choice of branches to formulate the TBA equations can be very 
important for its good definition and in particular for the future numerical applications. In this 
section we will fix a particular choice of branches in the kernels involved in the integral equations. 
This choice will be quite unique, with the following nice properties for the y-functions and the 
integration kernels: 

• They have only a minimal number of cuts, in general only a pair of cuts, which means that 
they obey an ordinary fusion procedure where all the intermediate constituents of a bound 
state but the first and the last cancel. 

• They are real functions of the spectral parameter u on the real axis. It fits well their physical 
meaning in TBA as of the ratios of densities of physical particles and holes. 

These properties will stem of course from the similar properties of integration kernels and free 
terms (with no convolutions) in the TBA equations (40)- (44). 

There are two natural possibilities to define x{u) compatible with (6). We define two functions 

They both solve (6). It is easy to check that with this choice of branches (7) reproduces the physical 
and mirror dispersion relations, correspondingly [ ]. They coincide above the real axes and have 
the following properties under complex conjugation 

;^ = xP^ , ^= l/x™'' . (48) 

Basically both representations (47) describe the same function, with the same Riemann surface but 
extended from the upper half plane to the plane with the cut {—2g,2g) for x^^, and to the plane 
with the infinite cut {—oo,—2g) U {2g,oo) for the function x™'''. One can say that they are two 
sections of the same Riemann surface. 

^Another possibility, consistent with the infinite length solution of [1], is to choose v £] — oo, —2g] U [2g, oo[. We will 
examine that possibility in detail in the next section. We thank G. Arutyunov and S. Frolov for the correspondence 
on this issue. 
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Figure 2: Structure of the cuts and conjugation paths on the mirror and physical sheets. 



We can plot them in Mathematica by running e.g. 
z=a+b I;xmr=l/2(z+I Sqrt [4-z"2] ) ;xpli=l/2(z+Sqrt [z-2] Sqrt [z+2] ) ; 
PlotSD [{Im [xph] , Im [xmr] +0.1},{a,-3,3},{b,-l,l}, PlotStyle->{Red , Yellow}] 

Notice that in the mirror ABA [14] (13) and (14) which we started from, the choice x™''' is 
employed [20, 14]. However, for the physical ABA of Beisert and Staudacher [ ] we only use the 
physical choice x^^. Thus to have a link with the ABA in the physical channel one should use the 
same definition (12) with 

Xj = x'^^iuj) , xf = x^^iuj ±i/2) , (49) 

in various free terms (with no convolutions) in the TBA equations. 

On the other hand, since all the kernels in the TBA equations are coming from the mirror 
theory, both arguments should be in mirror kinematics. Hence we specify in definitions (24) for the 
integration kernels the following branches'* 

, , x'^'Uu) - x'^'Uv) , l/x"'''(u) - x'^'^iv) , , 

r{u,v) = , . , b[u,v] = , . . (50) 

^x'^^iv) ^x'^^iv) 

With this choice of branches, it is easy to check that the kernels TZ^nm , ^nm entering our TBA 
integral equations (40)- (44) are all real! In the next section we show that the kernel involving 
the dressing factor, 2Snm, is real as a consequence of crossing, up to a simple square root factor 
which we identify there. Moreover, together with T^j'^^ ~ -^nm appearing in (44), it has very simple 
analytic properties. Namely, it has only four branch points for each of two variables, confirming 
the nice fusion property announced above. We will also present a simple integral representation for 
this combination. 



6.1 Reality and crossing 

One of the important consequenses of the crossing for the SU{2) x SU{2) principal chiral model 
considered in [i.i] was the reality of the function Yq corresponding to the single middle node in the 
finite size TBA equations. Here we show that exactly the same fenomenon is taking place in the 
present AdS/CFT TBA equations. 

Similarly to the Beisert-Eden-Staudacher physical dressing factor, the mirror S-matrix ought 
to be a pure phase. Let us here explain why this follows indeed in a simple way from the crossing 

^the same branches are used in [30, 31] 
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relation for the dressing factor. The same argument can be easily adapted to prove that the leading 
large volume y-functions found in [ ] are indeed real. 

We present schematically the mirror and physical sheets on the figures 2. They are naturally 
devided by cuts into three regions denoted by mi , m2 , ms and Pi,P2,P3, correspondingly. Since 
x^^{u) coincides with x'"^(u) in the upper half-plane the regions pi and mi are equivalent, pi = mi. 

Let us consider two points ua and va above the upper cut, i.e in the region pi = mi . Conjugation 
in the mirror sheet sends these points to ub = ua and vb = va (belonging to the physical sheet) 
while conjugation in the physical sheet maps them to uq = ua and vc = va (belonging to the 
mirror sheet). 

Notice that crossing condition relates the dressing factor with argument ub with the dressing 
factor at the point uc- More precisely, we have [-3] 

(7{uB,VB)cr{uc,VB) = —TTT ^ , x = x^'^iuB) , y = x^^{vB) (51) 

x^ — y l/x+ — 2/+ 

Notice also that we can now analytically continue both sides of this equality with respect to the vb 
root, in particular we can generate the crossing relation where is replaced by vc- Using again 
the (analytically continued) crossing relation to transform vb into vc we get 

o-{uc,vc) = a{uB,VB) , x = x'^^{ub) , y = x'^^{vB) (52) 

x+y 

Taking the complex conjugate of this expression and using the fact that the dressing factor is a 
pure phase on the physical sheet we get [14] 

{a{uA,VA)r = ^^^ r , x = x^''{ua) , y = x^''{vA) (53) 

X y+ a{uA,VA) 

Notice that we replaced x'^^{ua) and x^^{va) by their mirror counterparts because A is in the region 
pi = mi. Furthermore, in the left hand side, we explicitly wrote uc = ua and vc = va to recognize 
the explicit definition of the conjugated function on the mirror sheet. It is now clear that up to a 

simple factor of y'^f^^ the dressing factor in the mirror theory is indeed a pure phase function. 

More precisely the combination -^/^^^^(ti, v) is a pure phase in the real axis of the mirror sheet. 
The same kind of arguments can be used to prove the reality of the large L Y-functions of [i]. 



6.2 Integral representation 

We will show that the dressing phase on the mirror sheet admits some concise integral repre- 
sentation. Based on that representation we can explicitly see that it has very simple analytical 
properties. In particular, up to a simple multiplier, namely the simple square root factor identified 
in the previous section, we can clearly see that this dressing phase is indeed a pure phase function. 

6.2.1 A new representation of the dressing kernel in (mir,mir) kinematics 

We will start form the DHM integral representation [21] for a{x'^^{u it i/2),x^^{v it i/2)), 

a = exp [ix^^ + ix - - (54) 
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where x = =t V^i v ± i/2), 

I f dzi f dz2 1 1 . r(i?i;i - ^2 + 1) , , 

X{u,v) = - (p — (b -J— -J— log— : — - 55) 



kl| = l k2| = l 

and wi^2 = 5'(-2i,2 + l/-2i,2)- This representation is vahd for the physical kinematics and in particular 
for u,v in the region pi. Since pi = mi we can start with the same expression with x'^^{u) and 
x^^{v) replaced by x™'^(u) and x™''(f) for u and v in the region mi, above the upper cut. 

For the kernel S{u,v) = ^dvlog(j{u,v) appearing in our TBA equations we have 



S(u, v) 



29 r2g 
-2g J-2g 



(n^^^'^ {u, wi - iO) - (n, wi + iO)) g{wi - W2) x 
X (^7^(°^) iw2 - iO, v) - B^^^'^ {W2 + iO, v)) dwidw2 (56) 



where 



^ / N du , r(l — iu) / 1 \ 7 / „N 

g ^ -i^log = ^ i^,^ - _ + 1 . 57 

a=l ^ ' 

Let us briefly recall how to derive this representation for the dressing kernel from the integral 
representation (55). First the pole terms l/(2;i — x(n)) and 1/(2:2 — xiv)) are written as derivatives 
of log's which will give rise to the S's and 7^'s in this expression (in this section we often omit the 
lower indices of ;B's and 7^'s in which case they are equal to . . .11). The extra derivative to make 
a kernel out of the phase can also be transported to the log of gamma function by integration by 
parts and this generates the function Q. The integration contour around the unit circle in the 2:1^2 
variable is mapped to an integral from 1g to —2g slightly above the real axis and then back from 
—Ig to 2g slightly below the real axis for the variable w\^2- When w\ is above the real axis we 
have z\ = x^'^iwx) = x^'^'^{w\) but when we are below we have z\ = xP^(u;i) = l/x™"^(u)i). This 
explains why we get that combination of 7^'s and ^B's in the last formula. 

Now that we have transformed the original contour integrals into usual integrals in the real axis 
we can further replace the integration limits in this expression by =Foo because for \w\\ > 2g we 
have 7^^^*^) {u, wi — iO) — B^^^^ {u, wi +iO) = and similarly for W2- Hence we arrive at the following 
integral representation for the dressing kernel in the mirror kinematics 



00 roo 



S{u,v) = - / (n'^^^'>{u,wi-iO) -B^^^\u,wi+i0))g{wi-W2) X 
J —00 J —00 ^ ' 

X (n^^^'^ {w2 - iO, v) - {W2 + iO, v)^ dwidw2 . (58) 

Recall that this expression is derived for u and v in the region nii. The reason for which we cannot 
use this integral representation everywhere on the mirror sheet is the presence of poles of TZ under 
the integral at m = zizwi ziz i/2 — iO. To get rid of them we use the relations'5 7^(lO) = Ki - ^(^o) , 



^It is often useful to change from TZ^a to Z3's because the latter are much more regular than the former. In 
particular, since Im x'^" {u),lm > we can never have l/a;™"^(u) = x'^"{v) and thus B will be pole free when 

both variables are taken in the mirror sheet. Similarly, in the physical sheet, |a;P''(u)|, \x'^^{v)\ > 1 and again B is 
regular. Only when it and v are in different kinematics we should worry about regularity of the B functions. 



12 



7^(01) = Xi — ^gC^^) and then evaluate the integrals with Ki by poles using^ 



Ki{u - wi)g{wi -W2)B^°^'^ {w2 + i0,v)dwidw2 = -S^^^^ {u + i/2,v) - ^V^^'>{v) , 
J (u, wi + iO) g{wi - W2)Ki{w2 - v)dwidw2 = -B^^^^ {u, v + i/2) , 

(v) = dy log — . , . (59 

^ ^ 27r ^ - m/2) 

In this way, we get the following representation valid everywhere in mi , m2 , ms for both variables 

i 



where 



S{u, v) 



.^(11) (, 



oo fOO 



oo J — oo 



(n, wi — ^0) ) X 



Kg{wi - W2) f-8(°^) {W2 + iO, v) - (^2 - «0, t;) 



dwidw2 ■ 



(60) 



We see that the integrals can be combined to a contour integral around the cuts (— oo, —2g)\j{2g, oo)! 
This implies that we can write the result in a fashion similar to (55). Introducing 



ki|>i 



dzi 
27 



\Z2\>1 



dZ2 

27 



1 



{zi - x--(n)) {zi - x^'\u)) 



{Z2 - X^^'{V)) - X^^'{V)) 



log 



T{iui - iu2 + 1) 
T{iu2 - iui + 1) 



(61) 



with the integration going along the part of the real axes over (— oo,— 1) U (l,oo). Then for the 
physical dressing factor, analytically continued to the mirror in both variables we get the following 
representation 



1 - l/{x-y+^ 
1 - l/{x+y- 



cT™''"(u,i;) 



a{u, v) , a = exp [ix^+ + ix' 



(62) 



where x = x™''(ii), y = x'^"{v) and = x('W i */2, "U ± i/2). We see that the second factor a has 
the same properties under the fusion procedure on the mirror sheet as the physical dressing phase 
a had on the physical sheet - one simply replaces shifts by zbi/2 by =bin/2 for u and by zbim/2 for 
V in X- Note that X is a real function and thus it is a pure phase. Thus (7{u,v) is nothing but the 
dressing phase of the mirror theory! 

Finally let us present yet another interesting representation of the dressing phase in the mirror 
kinematics. It is easy to see that 7^(^°)(M,^^;) and n^°^\w,v), as functions of w, are regular below 
the real axis. Moreover B^^^\u,w) and B^^^\w,v) are regular on the whole complex plane except 
for the Zhukoswky cuts, see previous footnote. That implies that the terms with BB and 7^7^ 
in (58) vanish because for those terms we can deform the integration contour to +ioo and — ioo, 
correspondingly. For the remaining terms the integration with g can be done explicitly to yield 



2Snm{u, V) - 7e^^^) [U, V) + B^^^^} {u, v) 



-Kr. 



[U — V] 



a=l 



{u, w + ia/2) B^f^^ {w - ia/2, v) + B'J> {u, w - ia/2) B'Z' {w + ia/2, v) 



(10) 



j(01) 



(63) 



dw 



"in these formulae, w is the variable being integrated over in the last convolution 
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where we wrote the result already after fusion, i.e. for the dressing factor between magnon bound 
states n and m. Quite remarkably this combination of kernels, which is precisely the one appear- 
ing in the TBA equations contains no cuts apart from those at Im (u) = ibn/2 in the u plane 
and Im (v) = ±m/2 for the v variable, precisely as expected. This property was also noticed 
independently in [■>2]. 



6.2.2 A new representation of the dressing kernel in the (mir,ph) kinematics 



In this section we analyse the dressing kernel when a{u, v) when the first variable u takes values in 
the mirror sheet while the second variable v lives in the physical sheet. This is precisely the case 
for the free terms (without convolutions) in the TBA equations. For example, in (43) the term $ 
contains S{u) = Y\j (^{ujU^j) where U/j^j are the Bethe roots of the physical theory while u is in 
the mirror kinematics. The derivation of a nice integral re presentation for this dressing factor goes 
along the same lines as in the previous section. We find 



logS(u) 



(n, w + iO)- 7e(^°) (n, w - iO) 



log 



5(+)(n + iO) i?(+)(n-iO) 



Bi-){u + iO) 



log 



R(-){u-iO) 



(64) 



where R and B are defined like in (12) with x{u) = x™"'(u) and = x'P^{uj it As in the 

previous section this relation is derived in the region mi = pi and the next step is to transform this 
expression in such a way that it allows for a trivial analytical continuation to the full mirror sheet 
for the u variable. Actually the r.h.s. is not singular in mi,m2 (but not in rn-3), and thus should 
coincide with analytical continuation of the dressing factor. Next we recall that 



B'^+^Uj - i/2) = , R^-\uj + i/2) = 



(65) 



while log-B^ ^ and logi?*^"*"-* are regular in mi, 1712, m^. Assuming u to be real we can again drop 
BB and RR terms and convert 7^ to ;B as in the previous section to obtain 



logS(u) = B''^°\u,w + iO)*g*log 



R^+\u-iO) 
i?(-)(ti-iO) 



(66) 



^nnv ^ ^ , B^+\u + iO) ^ , B^+\u + iO) 

B'^^^'iu, w-iO)*g* log -^^^ — + Ki*g*log- ^ ' 



B(-\u + iO) 



5(-)(n + iO) 



The last term can be computed exphcitly, Ki*g* log = - log ^|-i(^+-/2) ~ ^ 

and in this way we obtain the following integral representation valid in the full mirror sheet mi, 
7712 and m3 



B(+)(M+i/2) 



log 5 



1 5^"^+ ^^-l, 4 



+ { B'^'^°\u,w + iO) *g *log 



(67) 



RM(u-iO) 



Ri-)( 



iO) 



+ ^(10)( 



u, w 



io) * g * log 



B(-\u + iO) 
B(+){u + iO) 
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Fusion is again trivial and yields 



m— 1 m — 1 r 1 



iu + i/2 + ik) ' ^2 ^.[-'"l ^^^^ 



+ (^g^ni^ (n, + .0) * g * log ^ , _ + Bl,,\u, w - ^0) * Q . log ^^^^ ^ , ^ j . 

Using the same arguments as in the previous section we could explicitly eliminate Q and one of the 
convolutions in this representation at the expense of introducing an extra infinite sum over a. In 
this way we could derive an alternative integral representation very similar to that in (63). 

In the next section, the reality property is discussed in further detail and in particular we explain 
why the y-functions which solve our integral equations are indeed real. 



6.3 Reality and Analyticity properties of y's 

We can easily check using the explicit large L solution for the y-functions presented in [1] together 
with the explicit representation of the dressing kernel derived in the previous section that all Y- 
functions are real when u is in the real axis^. To understand that this property actually holds for 
the y-functions even at finite L we should study the reality of several kernels in the TBA equations 
and also the reality of the free term (without convolutions) . If both are real then the exact finite 
L solution for y-functions will be also real. The reality property is of a particular interest for the 
future numerical applications of our equations which can be now done by iterations starting from 
the known large L solution. 

The most complicated equation to analyze is the one for the middle node, eq.(44), which contains 
the dressing factor in the (fused) kernel and in the free term. We will focus now only on this 
equation since the reality of all other equations can be checked trivially. Let us explicit in (44) only 
the ^^dangerous" terms: 

n-l 
2 

logy.„ = 25„^ * log(l + y,,J + J2 + ^A:) + 27^i2°) * log(l + y^j + . . . 

where the . . . stand for the rest of the terms, which are explicitly real. The reason why we also 
kept the last term as dangerous (i.e. potentially not real) will become clear below. 

Inside the kernel Snm the only non-real contribution comes from the square root of-ip("^) m 
(63) and the dangerous terms coming from the fusion of ^ are those in the first line of (68) so that 

^For the V-functions Yn and Y22 associated to the fermionic roots this property is true for \u\ < 2g. 
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we can re-write the dangerous terms in the r.h.s of (44) as^ 



logy.„ = -ipM*iog(i + y.^) + Viog^ (69) 



+ log 



B{-)(u-f) 



Q(n-iH±i)Q(n-i^ 



+ 27^i^°) * log(l + FaJ + . . . 



n+l\r^f„. I „■ 7 1—1 ' 
2 - 



Qiu + i^^)Q{u + i 



Now we notice that the first hne is nothing but the total corrected momentum (compare with the 
expression (2) for the corrected energy) which should vanish due to the string theory level matching 
constraint!'^ Thus, the only danger stems from both terms in the second line: in fact, the are not 
real (even though the kernel is real) but their combination will be shown to be real. 

The reason for the second term to be not real is that the function Yq^ contains pole divergencies 
on the real axis located precisely at the positions of the Bethe roots Uj. This can be seen from the 
free terms (containing no convolutions) in the TBA equation (42). We see that iff n = 2 then we 
do get singularities in the real axis coming from the zeros of R^~^{u + i/2) = which are precisely 
the Bethe roots uj. The zeros of this function induce, via this integral equation, the poles in Y[^^}^ 
Analyzing all other TBA equations in a similar way we can easily see that no other free terms give 
rise to poles in the real axis for any other ^-function. 

Now let us explain why the second line in (69) is explicitly real. The convolution in the presence 
of these poles should be understood as 2TZ^^2^ * log(l + Y[^^{u — iO)) which we can re-write as a 
principal value integral (which will be of course explicitly real) plus half of each residue of the 
singularities at the Bethe roots, i.e.^^ 

2<-) * log(l + Y^^iu - ^0)) = 27^i-) logd + J + log + ^^""'"^^^ + 



R(-){u-f)B(+){u 



m ^ 
2 ' 



The last term in this expression is not real. However it can be easily seen that it combines with 
the first term in the second line of (69) to give a real contribution! 

This concludes our check of reality of all the kernels and free terms in all TBA equations. The 
reality means that the y-functions solving these equations will be real, at least on some stretch of 
the real n-axis. 



7 Discussion and conclusions 

The integral equations we present are suitable for the numerical study. In the large L limit we can 
drop all convolutions containing the black nodes and recover in this way the large L solutions of 



To simplify the first line in (68) we use the identity 



B'--\u + i/2 + ik) 



BW(u + i/2 + ik) 



'^ B^+\u + i/2 + ik)R^-\u - i/2 + ik) 
- i/2 + ik)RM{u + i/2 + ik) 



= log 



B(-\u+f) 



B(-){u~f) 



Q{u + i!^)Q{u + i!^) 



^the gauge theory analogue of level matching is the ciclicity of the trace in the definition of local gauge invariant 
operators 

^''For n > 2 we also have poles for the corresponding Y-functions but they will lie away from the real axis. 
Notice that the residues at these singularities depend only on the prefactor of the log since when integrating by 
parts we get a log derivative which has always unit residue. 
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['] (we also checked this statement numerically). However, compared with the y-system equation 
in functional form these equations are of easy numerical implementation and the iteration from the 
large L solution to the finite L case is now accessible. This numerical approach is currently under 
investigation. 

In conclusion, we derived in this paper the system of non-linear integral equations of the TBA 
type describing, in principle, the spectrum of the states/operators in the full planar AdS/CFT 
system, including the low lying ones, such as Konishi operator. Not only these equations confirm 
our Y-system conjectured in [I] but they also give a practical way to the numerical calculation 
of the anomalous dimensions as functions of the coupling A. An alternative, usually numerically 
quite efficient, would be the derivation of the Destri-DeVega type equations along the guidelines 
presented in [12] for the SU{2) principal chiral field. In any case, a better understanding of the 
analytical structure of these equations is needed for the efficient numerics. 

A point which we do not completely understand in detail concerns the role of the so called fi- 
term contributions in the TBA equations. In particular we might need to pick extra contributions 
in (2) coming from further singularities which might arise in the l^^o functions. In the large L limit 
such extra terms could probably be identified with the Liischer's n term contributions. The role of 
these extra contributions, if they are present at all, needs to be further elucidated. 

One more unclear point concerns the underlying PSU{2,2\4:) symmetry of the problem. In our 
approach the starting point is the string theory in the light cone gauge where this symmetry is 
broken to SU{2\2)'^. It would be extremely interesting to understand how the full superconformal 
symmetry emerges in the TBA equations. 

Interesting questions yet to be considered concern the derivation of a full set of finite size 
Bethe equations for any type of excitations of the theory, again along the lines of [1 2] as well 
as the generalization of these TBA equations to another integrable example of the AdS/CFT 
correspondence, the ABJM duality [2G] (see [27, 28, 29] and references therein for the integrability 
related works on this theory). 

The set of TBA equations derived here should give us access to the full spectrum of AdS / CFT 
for any coupling. Hopefully it will help to understand deep physical reasons of the integrability of 
= 4 SYM theory. Knowing the exact results always helps understanding physics. 

Note Added 

After the work on this project was already finished the paper [.30] appeared where essentially the 
similar equations for the vacuum have been derived except the corner, fermionic nodes l2,±2- The 
corresponding equation 5.71 proposed in [30] appears to be incorrect. We derive here the correct 
equation and also propose the TBA equations for the excited states. 

Note Added for preprint arXiv:0902.4458v3 

When we were preparing a paper with the results of the new section 6, the second version of 
the paper arXiv:0904.4575 [■"'] appeared where a part of our new results, concerning the fusion 
properties for the mirror dressing factor, was established. We decided to update our old paper 
with these new results. We also stated more explicitly our preferable chose of the contours and 
branches in the integral equations. It is consistent with [20, 14] and agrees with that of [30, 31]. 
We also show that for that choice the y-functions for excited states have particulary nice analytic 
properties and are real. 

In the preprint arXiv:0902.3930v2 of ['JO], which appeared after our preprint arXiv:0902.4458vl of the current 
paper, the eq. 5.16 (equation 5.71 in arXiv:0902.3930vl), associated with ^2,2, was indeed recognized to be incorrect. 
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